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To construct an iV-representable time-dependent density-functional theory, a generalization to 
the time domain of the Levy-Lieb (LL) constrained search algorithm is required. That the action 
is only stationary in the Dirac-Frenkel variational principle eliminates the possibility of basing the 
search on the action itself. Instead, we use the norm of the partial functional derivative of the 
action in the Hilbert space of the wave functions in place of the energy of the LL search. The 
electron densities entering the formalism are iV-representable, and the resulting universal action 
functional has a unique stationary point in the density at that corresponding to the solution of 
the Schrodinger equation. The original Runge-Gross (RG) formulation is subsumed within the new 
formalism. Concerns in the literature about the meaning of the functional derivatives and the 
internal consistency of the RG formulation are allayed by clarifying the nature of the functional 
derivatives entering the formalism. 

PACS numbers: 



I. INTRODUCTION 

Density-functional theory (DFT) now provides the 
conceptual, theoretical, and computational framework 
for the study of the ground-state properties of a vast ar- 
ray of quantum-mechanical systems at all levels of aggre- 
gation from atomic to macroscopic. The foundations for 
the contemporary theory of chemical reactivity emerge 
naturally from DFT as well y . t| The essential ele- 
ments of DFT are the Hohenberg-Kohn (HK) theorems 
0, the Kohn-Sham (KS) equations ^j, the Levy-Lieb 
(LL) constrained search algorithm which together 

with the Harriman-Zumbach-Masche (HZM) construc- 
tion introduces A-representable densities into DFT; 
accurate approximate functionals 0; and powerful com- 
putational algorithms [To| . 

As defined through the LL algorithm, the density 
functional E[n] has a unique global minimum at the 
ground-state density within the space Af of all allow- 
able electron densities n(r). This variational prin- 
ciple of DFT stands in one-to-one correspondence with 
the Rayleigh-Ritz variational principle for the time- 
independent Schrodinger equation and provides the same 
generality to the derivation of the KS equations. 

Following the ground-breaking HK paper, a series 
of steps was taken towards a time-dependent density- 
functional theory (TDDFT) [U G3 which culminated 
in a more general formulation by Runge and Gross (RG) 
[l3l | . TDDFT is now being routinely applied to the calcu- 
lation of excitation energies of atoms and molecules , 
as well as various physical properties within the linear 
response regime and beyond it (see ref.0 for a 
survey of recent applications). 

In parallel to this success, discussions regarding the 



foundations of the theory continue to take place [l7| . 
In their original work, Runge and Gross |13| employed 
the quantum-mechanical action integral (hereon the RG 
action) to derive time-dependent Kohn-Sham e quat ions 
through the Dirac-Frenkel variational principle |l8j . It 
was later argued that the RG action led to para- 
doxes when calculating response functions because these 
must be causal, whereas second functional derivatives 
of the RG action were thought to be symmetric. This 
"symmetry-causality paradox" was resolved first by Ra- 
jagopal |2£|, who introduced an action based on the 
time path introduced by Jackiw and Kerman [2l|, and 
subsequently by van Leeuwen 22], who reformulated 
TDDFT replacing the RG action by a Keldysh action 
|23| . The RG, Jackiw-Kerman, and Keldysh actions are 
defined only for time-dependent w-representablc densities 
(TDVR) , and, regarded as functionals only of the density, 
are not stationary at the density of the solution of the 
time-dependent Schrodinger equation [l2l l22l 123. |25| . 
Such a lack of stationarity is a decided inconvenience, but 
even within TDVR TDDFT stationarity can be restored 
by recognizing that the density and external potential can 
be treated as independent functions [2(|. Nevertheless, 
for reasons analogous to those applying to the ground- 
state theory, it is important to generalize the definition 
of the action functional to hold for time-dependent N- 
representable densities (TDNR) . More explicitly, Mearns 
and Kohn have shown that small, time-dependent 
additions to the ground-state density need not be v- 
representable in first order. A suitable generalization 
can be effected by constructing a constrained-search al- 
gorithm for TDDFT analogous to the LL algorithm for 
DFT. 

Apart from restoring stationarity to the action in 
TDDFT, A-representability is important because, as in 
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DFT, accurate solution of the KS equations requires it- 
eration to self-consistency. The most convenient start- 
ing densities may well not be u-representable, nor may 
the densities be at intermediate stages of the computa- 
tional algorithms. It is then essential to have an action 
functional and KS potentials defined for iV-representablc 
densities both as a matter of principle and for practical 
reasons. 

In this paper we formulate an 7V-representable TDDFT 
based on the Dirac-Frenkel variational principle in which 
the RG action functional is stationary with respect to n 
at that unique n derivable from the solution of the time- 
dependent Schrodinger equation. We establish a one- 
to-one invertible map between all densities in a time- 
dependent generalization of Af and wave functions by 
use of the norm of the partial functional derivative 28] of 
the RG action in the Hilbert space of the wave functions. 
Insertion of that map into the action defines the action 
functional. The Runge-Gross formulation of TDVR is 
subsumed within this TDNR TDDFT, and the desired 
stationarity and generality are achieved. 

In Section II, we begin by reviewing two topics cen- 
tral to our later developments, the Dirac-Frenkel varia- 
tional principle and the action and its total and partial 
derivatives. Via Section II we introduce our notation for 
wave functions, operators, functional derivatives, Hilbert 
spaces, and more general function spaces. We also intro- 
duce the notion of mapping between abstract spaces as 
central to the formulation of TDDFT, following Dreizler 
and Gross for DFT 29] . In Section III, we recapitulate 
the RG formulation of w-representable TDDFT and show 
explicitly that its unnecessary limitation to the density 
generated by that v which enters the Hamiltonian de- 
stroys the stationarity of the action functional. We also 
provide an explicit explanation of why there are no in- 
consistencies in the functional derivatives entering the 
theory and why second functional derivatives of the RG 
action with respect to the density are not symmetric. Up 
to this point, our paper has concerned itself with the clar- 
ification of existing work on w-representable TDDFT. In 
Section IV, we turn to the problem of establishing a sat- 
isfactory 7V-representable TDDFT. We begin by stating 
a set of criteria that such a theory must meet. Next, we 
review existing proposals ( 26], 30]) and show that they 
do not meet all of the criteria. Finally, we develop the 
principal result of this paper, a constrained search algo- 
rithm which meets all of the criteria. We close with a 
brief summary of our results in Section V. 



II. BACKGROUND AND NOTATION 

The Dirac-Frenkel variational principle 

Consider a finite system of electrons and nuclei con- 
taining N electrons. Ignoring nuclear kinetic energy, 
keeping the nuclei fixed, and discarding the internuclear 
interaction energy as an irrelevant constant, the system 



Hamiltonian becomes 

H[v] =f + W + V[v] =H + V[v]. (1) 

In Eq.QJ, T is the electron kinetic-energy operator and 
W the electron-electron interaction operator. The opera- 
tor V[v] is the energy of interaction of the electrons with 
a time-dependent external potential v(r,t), 



V[v] 



dr v(r, i)n(r) 



(2) 



In Eq.|J2J, n(r) is the electron-density operator. v(r,t) is 
comprised of the potential energy of an electron in the 
fixed nuclear electrostatic potential plus that in a time- 
dependent potential generated by sources external to the 
system. For each time t in the interval (to,ii) under 
consideration, the r-dependence of v(r,t) must meet the 
conditions imposed by Lieb Q. In addition, we impose 
the requirement that 



«(r,t)-f0 , r-Tco , Vt e (to,h) 



(3) 



to eliminate irrelevant phase factors in the wave functions 
(see also ref . [3lJ ) . The time dependence of v must meet 
certain implicit integrability conditions discussed below. 
Such acceptable potentials lie in the space V. The space 
R 3 x (to,t\) is the support on which the elements v of 
V are defined. As indicated by our notation, V[v] is a 
linear functional of v, Eq.(|2J), and so, consequently, is 
H[v], Eq.Q. 

The wave-functions <I>(t) of the TV-electron system are 
time-dependent, normalized, antisymmetric functions of 
the N space and spin coordinates of the electrons, 



||$(t)|| = ($(«),$(*)) = 1 , Vte(t„,ti) 



(4) 



They satisfy the time-dependent Schrodinger equation 
(atomic units are used throughout), 



id t $(t) = H[v]$(t) 
Once the initial condition 

$(t ) = $o 
is imposed, &(t) is unique, 



*(t) = T L exp 



dt'H[v] 



(5) 



(6) 



(7) 



In Eq.||7|) Tr, is the time-ordering operator, later to the 
left. Eq.Q defines implicitly the conditions which v(t), 
<I>(i), and $0 must meet. In addition to those conditions 
which were specified by Lieb 0, &(t) must be differen- 
tiable in time. The set of such functions which are so- 
lutions of Eq.JSJ for all v in V form a Hilbert space 
They are supported in 3? on the space r, which is the 
product of (to,ti) with the configuration and spin space 
S of the iV-electrons, 



S x (t ,ti) 



(8) 
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All scalar products like the norm entering Eq.Q are de- 
fined on S. 

Eqs.© plus JHJ implicitly, and (0) explicitly, define a 
mapping Mi : V — > <&. Mi is surjective; $ contains no 
element which is not associated with an element of V [H| . 
That Mi is injective as well, i.e. one-to-one and therefore 
bijective or invertible, can be seen as follows. Suppose 
there is a v' and therefore a V' which yields the same <E> 
as solution of Eq.JSJ) as does v and V. Subtracting the 
two Schrodinger equations leads to 

V'-V^)$ = J dr[v'(r,t) - v(r,t)]n(r)<f>(t) = , 

which implies that v 1 must equal v under the conditions 
on the r-dependence of v required for the analogous proof 
for the time-independent problem (cf ref.|7j and p. 5 of 
ref. 29]). Thus Jl/f 1 exists and 3> <-> V is one-to-one. $ 
can then be regarded as a functional of v, $[v], or v one 
of $, u[<&]. 

Let us now expand the Hilbert space <fr to *ff which 
contains all functions which meet the conditions im- 
posed on $ including ^(io) = $0, except that the 
need not satisfy Eq.Q. The Dirac-Frenkel variational 
principle states that \& satisfies Eq.© if and only if 



id t - H[v] 



* = 



(10a) 



to define Frechet and Gateaux derivatives |28|-[32], which 
restrict the functions on which they are defined to a sin- 
gle Banach space (in the case of a Frechet derivative), or 
normed space (in the case of a Gateaux derivative). The 
Gateaux derivative has been regarded as a generalization 
of the concept of the partial derivative of a function [2^ . 
Similarly, the Frechet derivative has been regarded as a 
generalization of a total derivative |28j|. In our case, the 
properties of the action are such that taking derivatives 
only with respect to n meets the criteria for a Frechet 
derivative despite the fact that it is a partial functional 
derivative. In the following we shall use the terminology 
partial functional derivative to refer to derivatives with 
respect to a single function of functionals of more than 
one function. When, however, we map the potential in 
the action back to the density or vice versa so that the 
action becomes a functional only of a single function, we 
shall refer to the functional derivative taken with respect 
to that single function as a total functional derivative. 
The total differential could then be represented as a linear 
combination of partial functional derivatives times the 
corresponding differentials of the respective functions. 



III. V-REPRESENTABLE TDDFT 



1. v-representability and stationarity of the RG action 



(10b) 



The action and its total and partial functional derivatives 

We can now define the usual quantum-mechanical ac- 
tion functional A[^>, v] on the space $xV, 



A[V,v] = 







f 1 dt (% 


id t - H[v] 


/to v 





•) 



(11) 



v4[^>,w] is stationary only at $>[v] in iff with respect 
to variations d^ff, S^* taken at constant v, given that 
*(io) = $0 Vf £ *, and requiring as well that [34| 



(a*(ti),*(t x )) = 



(12) 



a less severe restriction than that of Eq. (jlObfl . The func- 
tional gradient of A[\&, v] along (v is fixed), 



G*. = dy,A[$,v] = id t -H[v] * , 



(13) 



thus vanishes in \& at yielding the time-dependent 
Schrodinger equation. Note the use in Ea. ((T3)) of dy* as a 
symbol for a partial functional derivative. We now clarify 
the nature of such a derivative. The action is a functional 
of two functions defined in two different spaces. Accord- 
ingly, it does not fit simple examples of functionals used 



We can recast the arguments of RG 13] as follows. Re- 
strict the argument ^ of A[*ff, v] in Ea. i|llfl to lie in <!?, the 
space of u-representable wave functions $[?/], defining an 
action functional, 



A[&, v}= I' dt ($, idt - H[v] $) 

J tn 



(14) 



on the space €> x V. Stationarity of A[iff, v] implies sta- 
tionarity of A[<5>, v] , i.e. that its partial functional deriva- 
tive vanishes, 



d$.A[$,v] = 



(15) 



since 3? C *Sf and the stationary point of A[^, v] is in 

A[<£, v] can be established as a functional of v alone by 
inserting in A[<5, v] that for which v' = v, 



A[v] = A[<$>[v],v] 



(16) 



The stationarity condition l|15|) then implies that the to- 
tal functional derivative of A[v] is — n, 



S v (v,t)A[v] = S v ( r tfA[$[v],v] = -n(r,t) 



(17) 



a generalization of the Hellmann-Feynman theorem |35j] . 
The total functional derivative (|17() does not vanish, ob- 
viously. It is only the partial functional derivative l|15|) 
which yields stationarity pq . 

To go on to the density-functional, A[n, v], requires 
establishing that the map M' 2 : $ — > J\f v , 



n(r,t) = ($(i),n(r)$(i)) 



(18) 
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is one-to-one and invertible. In Ea. (|18|) . n(r,t) is the 
time-dependent electron density, and the symbol Af v 
stands for the subset of all such v-representable densi- 
ties contained in A/", the time-dependent generalization 
of the space of densities of DFT All n(r, t) in Af and 
Af v obey the initial condition 



ra(r, t ) = n (r) = (<J> , n(r)$ ) 



(19) 



Eci. (|18fl defines what is meant by the phrase TDVR; a 
TDVR density is derivable via Ea. l|18|) from the solu- 
tion $ of the Schrodinger equation (jSJ) for some v in 
V. Demonstrating the invertibility of M' 2 directly, how- 
ever, is nontrivial. The HZM construction J8J shows that 
\I> — > Af is many to one. 

RG followed an alternative path. Substituting Eq.Q) 
into Eq.O defines a map M 3 = M\M 2 : V — > Af v . 
They then show by a pretty argument that M3 is one- 
to-one and invertible for all potentials v(t) which possess 
a Taylor expansion in time about to converging for all 
t £ (to,h). Van Leeuwen [25| has pointed out that it is 
sufficient for a Taylor series to exist about a set of points 
t{ G (to,ti) for which the radii of convergence overlap 
to cover (to,ti). M' 2 ~ x : Af v — > <fr can then be con- 
structed as M^Mi. Substitution of Mlf 1 , that is i>[n] 
into A[<f>, v], then yields the desired functional A[n, v}. 
A \n, y ] is stationary with respect to variation of n at fixed 
v 26], that is its partial functional derivative vanishes, 



d n A[n, v] = 



(20) 



It is important to recognize that A[n,u] is defined for all 
n generated via Eqs.Q and l|18l) from some v' , which can 
be varied independently of v. It is only at the stationary 
point that v' — v. 



2. Time- dependent Kohn-Sham equations 



Substitution of both M' 2 



-1 and M 3 1 , i.e. <&[n] and 



v[n], into A[$, v], then yields a functional A[n] of n only 
(for a given initial state |3(|), the RG action functional. 
One thus has the option of using n or v as the indepen- 
dent variable in the functional. Van Leeuwen [2f| gives a 
simple and elegant argument for the construction of the 
TDKS equations from A[n] without invoking stationarity 
in Af v . Switching now to A[n] from A[v], we carry out a 
Legendre transformation to 

B[n] = A[n] + [ dt f dr v(r, t; [n])n(r, t) . (21) 

From Eq.(J2IJ, it follows that 

S n (r,t)B[n]=v{r,t). (22) 

The TDKS equations [25| follow from J^J. 

Consider a system of non-interacting electrons de- 
noted by subscript s which move in an external potential 



v s (r,t), starting from a single determinantal state $os 
at to- v s is a functional of their electron density, w s [n s ]. 
The HZM construction Q allows identification of n s (r, t) 
with the density of the interacting system, 

n s (r,i) =n(r,t) , Vr,f e (t ,h) . (23) 

Thus v s can be regarded as a functional of n. Combining 

5 n (T,t)B s [n] = v s (r,t) (24) 
with Eq.J221) leads to 

v s (r,t)=v(r,t)-S n(r , t) (A-A s ) . (25) 
The usual rearrangements in A — A s in turn lead to 

V s (r,t) =v(r,t)+v a (r,t)+v xc (r,t) , (26) 



v H (r,t) = J dr 



, n{v',t) 
|r - r'l 



u X c(r,t) = -(5 n ( r ,()4xcW 



(27) 



(28) 



A xc [n]= dt{[(y,id t y)-(V s ,id t y s )} (29) 

Jt 

-[(T-T S )-(W-W H )]} . (30) 

T is the kinetic energy and W the energy of electron- 
electron interaction of the interacting electrons in state 
$ [n] . T s is the kinetic energy of the noninteracting elec- 
trons in state <!> s [n] . W u is the Hartree approximation to 
W using or equivalently $ s [n]. 

It is at this point that concern about the meaning of 
the functional derivative defining u xc , Eq. l|28[) . arises in 
the literature [30 |H HI ^ . Since in Section IV we 
shall base our development of 7V-representable TDDFT 
on the RG action and since the above-mentioned concern 
raises doubts about the validity of doing this, we now 
summarize the debate and show why the RG action is 
perfectly suitable for the developments of Section IV. 



3. The symmetry-causality dilemma 

Taking the functional derivative of Ea. (|26|) with re- 
spect to n results in |37j 

X - 1 {v 1 t ] v',t')= X - 1 {v,t-v' 1 t')+f{v,t-v',t') , (31) 

where 



, Sn(r, t) 

X(r,t;r',t')- 



Sv{r',t') ' 



(32) 
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Xs{r,t;r',t') = 



5n(r,t) 
'6v s (r',t') 



(33) 



, ,s_ S[v H + v xc ](r,t) _ 5{t-t') i <5wxc(r,i) 
j ^r, t, r , r ) — 



6n(r',t') 



r — r 



Sn(r',t') 

(34) 

From Eas. (|18l) and Q), the well known retarded char- 
acter of the time dependence of the susceptibilities x 
and Xs follows; they vanish if t' > t. Their inverses 
X _1 (r, t; r', t') and x7 1 ( r , ^ r', *') entering Ea. l|3"T|) are re- 
tarded as well. Yet 



/xc(r,t;r',i') 



Sv xc (r,t) _ S 2 A xc [n] 
Sn(r',t') ~ 5n(r',t')6n(r,t) 



(35) 



is formally a second derivative of ^4 X c[ n ] according to 
Ea. (|25|l . Van Leeuwen [52, assumes that, as / xc is 
a second functional derivative, it must be symmetric in 
r, t and r', t' . Such symmetry is inconsistent with the re- 
tarded nature of x -1 and x^ 1 in Ea. lj3"T|) . Van Leeuwen 
[22I EH describes this inconsistency as a "paradox" and 
develops TDVR TDDFT from the Keldysh action in- 
stead of the RG action to avoid it. The second functional 
derivatives remain symmetric on the Keldysh time con- 
tour but become retarded when map ped into real time. 

Gross, Dobson, and Petersilka on the other hand, 
suppose that Ea. lj3Tl) holds and that / xc (r, i; r', £') must 
be retarded and not symmetric in r,i and r',t'. T hey 
then conclude by supposing from Schwarz's lemma [33 
that (1) / X c(r, t; r', t 1 ) cannot be a second functional 
derivative and that (2) the exact v xc [n] cannot therefore 
be a functional derivative. They conclude further that 
this in turn is in contradiction to the principle of station- 
ary action which leads to functional derivative. 

To complicate matters further, Harbola and Baner- 
jee |38| have argued that there is no symmetry-causality 
dilemma because while ~x i s causal, x _1 is symmetric. 
Amusia and Shaginyan |41|, while not disagreeing with 
this conclusion, have argued that, in contrast, it is pos- 
sible to construct a causal y _1 as well. Harbola ^2 nas 
responded that reference [4l|] itself implies a causality in 
the potential as a functional of the density, van Leeuwen, 
however, has argued that x must be rigorously causal 
in analogy with the properties of discrete lower triangu- 
lar matrices |25j . an argument which does not take into 
account the fact that x _1 is n °t a smooth function of 
t — t' but contains both a delta function and the sec- 
ond derivative of a delta function at t = t l+ . We show 
in Appendix A that x _1 (^ — O is causal, consisting of 
those singular functions at t = t' + plus a smooth causal 
function of (t — t'), so that the dilemma remains. 



4- A way out of the dilemma 

We conclude that in the context of TDDFT at the 
RG level, x _1 is causal. The most forceful argument 



that the causality of x" 1 imposes a symmetry-causality 
dilemma via Ea. (|31l) is that of Gross, Dobson, and Peter- 
silka [l9|. The flaw in their reasoning is the supposition 
that Schwarz's lemma can be applied to the functionals 
of TDDFT. Throughout all of density- functional theory, 
the Frechet definition j2^| of the functional derivative was 
implicitly used. In the present instance, the functional 
derivatives of Eas.l|2n jl .p2 )l .ipi )l and are all Frechet 
derivatives. Taking a second derivative simply involves a 
single iteration of the Frechet operation [28| . For the first 
derivative to exist, both the functional and the function 
space must meet smoothness criteria. The first deriva- 
tive remains a functional, which for the second derivative 
to exist, must remain smooth. This condition is implic- 
itly assumed for v xc [n] in all of DFT and TDDFT, and 
we presume it here as well. We conclude that all of the 
second functional derivatives encountered in TDDFT are 
perfectly well defined iterations of the Frechet derivative 
operation. These include 



XM; r',t') 



5 2 A[v] 



6v(ji,t')6v(r,t) 



(36) 



and Xs as well as / xc - All have a retarded dependence 
on t and t' and are decidedly not symmetric in r, t and 
r', t' . Similarly, x _1 and x^ 1 can be expressed as second 
derivatives, e.g. 



X - 1 (r,<;r',0 = 



8 2 B[n] 



6v(r,t) 

Sn(r',t') ~ Sn(r',t')Sn(r,t) 



, (37) 



have retarded time dependence (Appendix A), and are 
not symmetric in r, t and r', t'. 

We therefore agree with the main conlcusion of Amu- 
sia and Shaginyan an d Harbola and Banerjee 
S3 that there is no conflict between the symmetry 
and the causality. However, the way out of the dilemma 
is not by finding symmetry in the inverse response func- 
tions, but by recognizing that second-functional deriva- 
tives need not be symmetric functions of the time vari- 
ables. To understand how this asymmetry can come 
about in a second functional derivative, consider that 
functionals are defined on three levels. First, there is the 
space on which the functions are defined; second, there is 
the function space on which the functionals are defined; 
and third, there is the definition of the functional. For 
example, A[v] is defined through Ea. (|ll|l and the map 
M\ , on the function space V within which the potentials 
v are supported on K 3 x (to, ii)- For the total functional 
derivative 6A[v]/Sv(r,t) to exist and equal — n(r, t), first 
V must be smooth enough that variations Sv(r,t) exist 
which can be taken continuously to zero. Following Lieb 
, we have defined V for this to be the case. Second, the 
functional A[v] must be smooth enough that the resulting 
variation in it, cL4[v], exists, is linear in Sv(r,t), and goes 
continuously to zero with Sv(r,t). A[v] meets that crite- 
rion. The functional derivative 6A[v]/6v(r,t) is then de- 
fined through Frechet 's theory of linear functionals [28| . 
Similarly, for Sn(r, t)/Sv(r', t') = 5 2 A[v}/6v(r' ,t')Sv(r,t) 
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to exist, n(r,i) need only meet the smoothness criterion 
as a functional of v, which it does through the definition 
of Af v . 

The requirement for the applicability of Schwarz's 
lemma, that the second derivative be invariant with re- 
spect to interchange of the order of differentiation, is that 
the first level of support, the space on which the func- 
tion is defined, be unchanged by the first functional dif- 
ferentiation. That is not the case here, and Schwarz's 
lemma docs not apply. In A[v], v is supported on 
K 3 x (<o,£i), but in n[v], the first derivative, v is sup- 
ported on R 3 x (to,t) precluding the applicability of 
Schwarz's lemma (see also the discussion in Appendix 
B). lit' > t in 5 2 A[v]/Sv(r' ,t')Sv(r,t), it must vanish, de- 
stroying symmetry while remaining a well-defined second 
functional derivative |2q|. In the Keldysh action func- 
tional used by van Leeuwen [22I l25j . the time-ordered 
contour on which the action is defined provides the sup- 
port for the time-dependence of the potential v. Func- 
tional differentiation of the Keldysh action does not mod- 
ify this support, and so the second functional derivative 
remains symmetric in that support. Transformation from 
the Keldysh time contour back to real time introduces the 
asymmetry without changing the fact that a second func- 
tional derivative was taken. Thus, van Leeuwen has, in 
effect, proved that second functional derivatives need not 
be symmetric. We conclude that all functional deriva- 
tives in the RG formulation are well defined, both first 
and second, and that Ea. (|31[l . a relation among second 
functional derivatives, contains no inconsistencies. One 
thus has a choice - one can base TDVR TDDFT on the 
RG action or on the Keldysh action. How to general- 
ize the Keldysh action so as to provide a basis for TDNR 
TDDFT is not now clear. Accordingly, we choose to base 
our development of TDNR TDDFT on the RG action. 



Previous work 

Apart from formulations applicable to special classes 
of potentials |43j , there are two proposals for the formula- 
tion of NR TDDFT. That of Kohl and Dreizler [U does 
not meet criterion 2.) and, as a consequence, cannot meet 
criterion 4.) as well. That of Ghosh and Dhara |2g does 
not produce A^-representablility, only w-representability. 
Their Theorem 4 can be restated as defining the map 
M GD : Jsf -» 

= ARG{STAT^^ n B[^}} . (39) 

However, since in Ea. l|39|l one searches only for a sta- 
tionary point of B[^f] in one is allowed to relax the 
subsidiary condition l|38|l by a Legendre transformation. 
Ea. (|39|) then becomes 

*[n] = ARG | STAT 3, - dt J dr A(r, t)n(r, t) 

Now, the Lagrange multiplier A(r, t) will exist if and only 
if n(r,t) is u-representable, in which case Eq. l|4U|l be- 
comes 

¥[n] = ARG {STAT* A]} (41) 

with A C V, a potential. Thus, the map defined by 
Ea. (|39|) is identical to the map M'f 1 :Af v ^& defined 
by RG. Eas. (|3"9"jl - l{4*Tl) should therefore be rewritten with 
$[n] replacing 9[n]. What Ghosh and Dhara have ac- 
tually accomplished is to find a simpler and more direct 
proof of w-representability than the original proof of RG. 

Our approach 



IV. iV-REPRESENTABLE TDDFT 



The HZM construction |8( establishes that at each time 
t, there is an infinite set of wave functions &(t) which 
yield any preset n(r, t) in M via the mapping : *ff — > 



n(r,i) = (*(*), fi(r)*(t)) 



(38) 



with n(r, to) = ?to( r )> Ea. (|I9|l . The task in constructing 
an A^-representable TDDFT is to select a single member 
of that set so that becomes one-to-one and invertible, 
i.e. to find M 2 _1 : M -> M^ 1 should meet the follow- 
ing four criteria. I.) It should be universal; 2.) it should 
require searching only in *S> and not in '3/ and V; 3.) 
it should subsume the mapping M!^ 1 of v-representable 
TDDFT; and 4.) it should provide a stationarity princi- 
ple. 



We note that the stationary point $[n] of B[^] in 
Ea. (|39|) is unique in the subspace <& n — > n) of \& 
for n C Af v . The partial functional derivative of A[fy, v], 
its gradient in vanishes uniquely there, 



idt-H 



d**A[*,v]) Vtn = d 9 .B[9]) n = 

= ; ncA/\,,tf = $[ra] . (42) 

Thus the magnitude squared of the gradient, 



dt(dy*A,d<s,.A) v , 



dt{d^,B,d**B) n (43a) 



has a unique minimum there as well. As the search can 
be restricted to normalized ^'s without penalty, it follows 
from eo. (|43a|l that 



r dt{dy.A,d^A) vn = r dt(*,[id t -H]H 

Jto Jto V 

( 43b ) 

holds as well because of the consequent hermiticity of idt- 
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On the other hand, no such minimum can exist in the 
magnitude of the gradient for an iV-representable n C M 
which is not v-representable. A similar situation exists in 
time-independent DFT. A minimum exists in the func- 
tional E[^\ = ($!,H$n for * — » n if and only if n is 
v-representable. If we suppose that a minimum exists 
under the constraint of fixed n for n not u-representable, 
the constant can be eliminated by a Legendre transfor- 
mation. The Lagrange multiplier then simply adds to the 
external potential contradicting the hypothesis that n is 
not w-representable as in the arguments associated with 
Ecis. H39fl - i|41fl . For a general ./V-representable n, there is 
only an infimum in both (^/,H 1> i>) and (^>,H^) at the 
same point. The Levy-Lieb constrained search algorithm 
makes use of this infimum to define the density functional 
for 7V-representable densities: 



E[n] = INF^^ n E[^} ; 



(44a) 



V. SUMMARY 



An iV-representable time-dependent DFT has been es- 
tablished, and a time-dependent analog of the Levy-Lieb 
constrained search algorithm has been proposed. The 
central quantity in this search is the norm of the par- 
tial functional derivative of the Runge-Gross action in 
the Hilbert space of wavefunctions. The proposed con- 
strained search meets all of the requirements we pose: it 
is universal, requires searching only in one Hilbert space, 
subsumes Runge-Gross v-representability, and provides a 
stationarity principle. 



APPENDIX A: CAUSALITY OF 1 AND x _1 



*[n] = ARG {INFy-mEl®]} 



(44b) 



An analogous constrained search algorithm can be con- 
structed for TDDFT from the magnitude of the gradient 
0: 

9[n] = ARG^INF^^ n J* dt{d^A,d^A) v Jj 

= ARG ^INF 9 ^ n dt (9*„B, 9*.-B) n | (45a) 
= ARG {lNF^^ n 1 dt [id t - H] 2 V 



A[n, v] 



to 



dt (v[n], [id t 



(45b) 



We note that the quantity over which the search in 
Eq. (|43al) is done corresponds to the time- integral of the 
McLachlan functional for H ^jjj widely used in formu- 
lations of semiclassical dynamics ^(J. The proposed 
constrained search algorithm expressed in Eas. l|45al45b|) 
meets all of the criteria imposed above: 1.) It is univer- 
sal, not involving v. 2.) It requires searching only in Vt. 
3.) It subsumes w-representable n for which the infimum 
becomes a minimum and yields the condition 



(idt-A)' 



<J = s.t. <f 



(46) 



which yields the same \fr as Ea. H39|) or, ultimately Eq.©. 
Finally, 4.) it provides a stationarity principle since n 
in l|45bfl can be varied independently of v, and the 
corresponding partial functional derivative vanishes via 
Eq.fl3, 



d n A[n, v}=0 



(47) 



Stationarity in TDDFT plays the role that minimality 
does in DFT regarding error reduction. 



As stated in Section III. 3 and III. 4, there is a substan- 
tial spread of opinion in the literature with regard to the 
time dependence of xj 1 and x~\ differing as to whether 
it is causal or symmetric. We argue here that it is un- 
equivocally causal, with local singularities at t — t' + . 
It is easiest to see this explicitly for the Xs of the uni- 
form electron gas, which has the form Xs(|r — r'|, t — t'), 
from space-time uniformity. Accordingly, it is diagonal- 
ized by Fourier transforming on space and time yielding 
the eigenvalues Xs( c li w ~ i^)> 5 I 0. The wave- vector 
q is introduced by the Fourier transform on r — r', the 
frequency lu is introduced by that on t — t' , and S is in- 
troduced by the causality of Xs, Xs(\ r — t — t') = 0, 
t > t'+. 

The explicit form of Xs(q, w — iS) is known ^t|- When 
continued to the entire complex angular frequency plane 
in the process of inverting the Fourier transform on time, 
its only singularities are a second-order pole at infinity 
and bounded branch cuts just above the real axis. For 
q < 2k F (kp is the Fermi wavenumber), there is one 
bounded branch cut at z = ui + iS, with |a>| < -^(2kpq + 
q 2 ); for q > 2k F there are two, with -^{~2k F q + q 2 ) < 



M < ^k F q- 
the branch cuts. 



Xs (q, has no zeros away from 



Xj 1 is also uniform in space and time and therefore 
diagonalized by Fourier transformation. It's eigenvalues 
are simply 



(Ai) 



The second-order pole in \s (q, at to = oo yields the 
following behavior in xj 1 (q, w) at oo, 



x s ^q,^) 



A s (qV 2 +B s (q) + ^M + ... (A2) 



\U)\ — > oo 

6^0 



8 



B s (q) 

CM 



m 
nq 2 
2E F 
5 n 
IQEpq 2 
Ylhnm 



5h 2 q 2 



8mE F 



3 + 35 



h 2 q 2 
8mE F 



where n is the number of electrons per unit volume and 
m is the electron mass. Thus x^ 1 has the form 

Xj'i^t- t') = -A s (q)S"((t - t')+) + B s (q)5((t - t')+) 

+x; 1 (q,t-i')'(A3) 

In Eq. (|A3|) , 6" is the second derivative of the delta func- 
tion. (xj 1 )' arises from the branch cut(s) and is rigor- 
ously causal because the locations of the branch cuts in 
X~ 1 (q, z) are identical to those of Xsil, z ), being in the 
upper-half z-plane. 

The principal change in passing from Xsf^^) to 
x(q, of) for the uniform electron gas is that the free- 
particle excitations are replaced by quasi-particle exci- 
tations which have finite lifetime except at q = 0. This 
causes the branch cuts to extend to infinity, but causes 
no change in the formal structure of x -1 (qi t — t') which 
is given by i|A3(l with modification of A s (q) to A(q), etc. 

For a non-uniform extended system for which the ex- 
citation spectrum forms continua, be the system ordered 
or disordered, there is no change in formal structure of 
X7 1 (r, r';t — t') and x _1 (r,r';t — t'). Each contains 
the local contributions S"((t - t')+) and 5((t - t') + ) 
as well as non-local retarded contributions. For finite 
systems which have at least one discrete excitation as- 
sociated with a transition from the ground state to a 
bound excited state, there is a change. Each eigenvalue 
of Xs( r i r '; w ) or x( r i r '\ u ') switches from +oo to — oo as 
the pole at z = to + iS with huj equal to that discrete 
excitation energy is crossed. This forces the existence 
of a zero between discrete excitation energies or between 
the highest discrete excitation energy and the continuum 
threshold. Each such zero gives rise to a pole in the cor- 
responding eigenvalue of x _1 ( r > r '> ° r X7 l ( T > r '> z ) at 
the same z. Upon Fourier transform to the time-domain 
X _1 (r, r'; t—t') and x7 1 ( r : r '; t—t') each contains a causal 
contribution from the pole which oscillates with angu- 
lar frequency corresponding to the excitation energy for 
t > t' and vanishes for t' > t. 

In conclusion, the causality of x( r 7 r ';^ — an d 
Xs(r, r';t — t') forces the eigenvalues of x( r : r '; w ) an d 
Xs(r, r';o;) to have singularities only in the upper- half 
complex-frequency plane. The corresponding eigenval- 
ues of x -1 ( r ; r ';w) and XJ l { r i r 'j w ) can therefore also 
have singularities only in the upper- half plane apart from 
the second-order pole at oo. This arises from the fact 
that u> enters x an d Xs only in the combination u> — iS, 
6 I 0, which does not change when their eigenvalues are 
inverted to obtain x _1 a nd xj 1 - Those quantities must 
therefore always be of the form 

X\t-t') = A6"((t-t')+)+B6((t-t') + )+x- 1 (t-t'y , 

(A4) 



where x 1 {t — f)' is nonlocal and causal in time. 



APPENDIX B: SECOND-FUNCTIONAL 
DERIVATIVE ASYMMETRY 

In section III below Ea. H37|) . we have pointed to the 
modification of the support of u(r, t) in n[v] by the first 
functional derivative of A[v] with respect to v as the ori- 
gin of the asymmetry of its second derivative, x- Al- 
ternatively, one can preserve the support on which the 
function v(r,i) is defined, but then the second level of 
definition, the function space on which the functional is 
defined, must change. Consider, for example, the map 
Mi : V -> *, Eq.Q), which, together with Eq.||T8J| de- 
fines the map M3 : V — » Af v , n = n[v}. A more explicit 
expression of that map would be 



n = n[H[v\] 



(Bl) 



according to Eq.JJJ, in which v(r, t) is supported on 
(to,t). However, Eq.Q can be rewritten as 



$(t) = T L exp 



-i / dt"H t [v] 

J t 



(B2) 



where 



H t [v] = H[v] , t' S (t ,t); 
= , t' G t,t x ) 



(B3) 



Thus, by changing the operator space on which the argu- 
ment of the functional, now the operator TL t [v\ 1 is defined, 
we have formally restored the support of v to (£o,ti). 
However, that does not eliminate the asymmetry; it triv- 
ially shifts the location of its origin, viz 



5n(r,t) _ 6n(r,t) 6H t 



Sv(r',t') SH t Sv(r',t') 



= 



t' > t 



(B4) 
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